## Import glacier data
GRTE_glaciers2019 <- st_read('data/Glaciers_2019.shp')
## Reading layer `Glaciers_2019' from data source
## `C:\Users\wreis\OneDrive - Colostate\MS Research\R_Git\TetonProject\TetonProject\data\Glaciers_2019.shp'
## using driver `ESRI Shapefile'
## Simple feature collection with 12 features and 3 fields (with 1 geometry empty)
## Geometry type: MULTIPOLYGON
## Dimension: XY
## Bounding box: xmin: 511960.3 ymin: 4841322 xmax: 519221.9 ymax: 4860846
## Projected CRS: NAD83 / UTM zone 12N
outlines <- mapview(GRTE_glaciers2019)
norm_rad_2019 <- raster('data/normSolarRad_2019.tif')
raster <- mapview(norm_rad_2019)
## Warning in rasterCheckSize(x, maxpixels = maxpixels): maximum number of pixels for Raster* viewing is 5e+05 ;
## the supplied Raster* has 1867764
## ... decreasing Raster* resolution to 5e+05 pixels
## to view full resolution set 'maxpixels = 1867764 '
## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded ellps WGS 84 in Proj4 definition: +proj=merc +a=6378137
## +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null
## +wktext +no_defs +type=crs
## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded datum World Geodetic System 1984 in Proj4 definition
outlines+raster